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Establishing that a signal of new physics is undoubtly supersymmetric requires not only the 
discovery of the supersymmetric partners but also probing their spins and couplings. We show that 
the sbottom spin can be probed at the CERN Large Hadron Collider using only angular correlations 
in the reaction pp — * bb* — » bb$T, which allow us to distinguish a universal extra dimensional 
interpretation with a fermionic heavy bottom quark from supersymmetry with a bosonic bottom 
squark. We demonstrate that this channel provides a clear indication of the sbottom spin provided 
the sbottom production rate and branching ratio into b\\ are sufficiently large to have a clear signal 
■ above Standard Model backgrounds. 
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INTRODUCTION 



Supersymmetric models are promising candidates for physics beyond the Standard Model (SM), despite the present 
lack of direct experimental evidence on supersymmetry (SUSY). The CERN Large Hadron Collider (LHC) has a large 
i-^ ' reach for the discovery of SUSY [l[ that largely relies on the production and decay of strongly interacting new particles 
i.e. squarks and gluinos @, S|- Establishing that a signal of new physics at the LHC is indeed supersymmetric requires 
O !* not only the discovery of the new supersymmetric partners but also probing their interactions and spins [J, [5|, [y, [7[ • 

Previously, the squark spin has been studied through the long decay chain q — > x° ~> & ~~ * X? which is also used to 
measure its mass [3, E| or in conjunction with the gluino spin analysis in the decay chain g ^ b ^ X2 ^ & ^ Xi 0- 
In this work we probe the potential of the CERN Large Hadron Collider (LHC) for unravelling the bottom squark 
^ ■ spin using angular correlations in the short decay chain pp — > bb* — ► bb-ft?, analogously to what has been done 
for the analysis of the slepton spin [1, A nice feature of this reaction is that only sbottoms and the lightest 
SUSY particle (LSP) takes place in it, providing a further check of sbottom spin obtained in the long decay chain 
studies. Unfortunately, this analysis can not be extended straightforwardly to light flavor squarks due to large QCD 
, backgrounds. 

■ To determine the spin nature of the sbottoms at the LHC we compare the SUSY sbottom production and decay 
I chain with another scenario where the new intermediate states have the same spin as the SM particles and leads to 
the same final state bb^r- Such a model are Universal Extra Dimensions (UED) Q where each SM particle has a 
heavy Kaluza-Klein (KK) partner which can mimic the SUSY cascade decay, provided we employ the mass spectra 
extracted from the decay kinematics match Q . Here we are not focusing on UED searches but we use UED only for 
comparison with the SUSY predictions. 

N 

There are many observables which we can use to discriminate 'typical' UED and SUSY models, like the production 
rate or the mass spectrum. Nevertheless, at the LHC we measure only production cross sections times branching 
ratios, and the UED as well as the SUSY mass spectra are unlikely to be what we currently consider 'typical'. On 
the other hand, spin information is generally extracted from angular correlations. Therefore, we base our analysis 
exclusively on distributions of the outgoing SM b quarks as predicted by UED and by SUSY. We demonstrate that this 
final state bb^T provides a clear indication of the sbottom spin provided the sbottom production rate and branching 
ratio into 6x? are sufficiently large to have a clear signal above SM backgrounds. 



'Electronic address: 
T Electronic address: 



aalvesOfma. if. usp.br 
eboli@fma .if.usp.br| 



2 



II. UED INTERACTIONS AND PARAMETERS 

We assumed one extra dimension with size R, where all SM fields propagate [1, [9[ , leading to a tower of discrete KK 
excitations for each of the SM fields (n — 0). In this scenario, the 5-dimensional wave functions for an SU(2)-doublet 
fermion are 



71=1 



On the other hand, the roles of the left and right handed n-th KK excitations are reversed for SU(2) singlets. Just 
like in the MSSM, the spinors of the singlet (q) and doublet (Q) KK-fermion mass eigenstates can be expressed in 
terms of the SU(2) doublet and singlet fields ipd,s 

= cosa (n Vd n) + sino»4 n) , 
g< n > = sina (n) 7 5 ^' l) - cos a (n) 7 5 ^ n) ■ (2) 

In general, the mixing ang le qW is suppressed by the SM fermion mass over the KK-excitation mass plus one-loop 
corrections except for the top quark due to its large mass. 

tan2a< n > = ^ ^ (3) 

The non-degenerate KK-mass terms 6m^ n ' contain tree level and loop contributions to the KK masses, including 
possibly large contributions from non-universal boundary conditions. 

The neutral KK gauge fields will play the role of neutralinos in the alternative description of the sbottom production 
and decay. Just as in the SM, there is a KK-weak mixing angle which for each n rotates the interaction eigenstates 
into mass eigenstates 

7^ = cos 4") + sin WQ , 

4") = - sin W fl W + C o S W . (4) 

The n-th KK weak mixing angle is again mass suppressed 

tan 2* W = __ v l99rll (5) 

(5m$ 3 )*-(SmP)*+v* (g 2 -g 2 Y )/A 

where 8m^ contains tree level as well as loop corrections to the KK gauge boson masses. Generally (Sm^) 2 — 
(dm^) 2 3> v 2 (g 2 — gy) [§] and the lightest KK partner is the with basically no admixture from the heavy 

We only considered the first set of KK excitations to formulate an alternative interpretation of production and decay 
of sbottoms at the LHC, using the UED decay — > bj^ to mimic a sbottom decay. The relevant interactions for 
this decay are 



g(°V (Y d cos aWp L + Y s sma^P R ) Q (1) - <z ( °V (V d sin a w P L + Y s cos q (1) 7W . (6) 



In general, the KK partners of the SM particles do not have a mass spectrum similar to what we expect in SUSY, 
however, we imposed that the first KK excitations have the same mass as the SUSY particles. Assigning the LSP 
mass to the Lightest KK Particle (LKP) and the Next-LSP mass to the Next-LKP mass fixes the KK-weak mixing 
angle by means of Eq. ([5]) 

0« = l - arctan ( - 2 2 ^ . (7) 

2 \^m 2 NLSP -m 2 LSP )J 
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test point 


particle 


mass (GeV) 


branching ratio 


S5 


Xi 


122. 


stable 


S5 


bi 


633. 


Br(6j bxi) = 26.6% 


S5 


b 2 


663. 


Br(6 2 -> 6x?) = 78.8% 


SPSla 


xi 


97. 


stable 


SPSla 


bi 


517. 


Br(6i -> bx°) = 4.4% 


SPSla 


b 2 


547. 


Br(6 2 -> 6x?) = 29% 


LI 


xi 


97.8 


stable 


LI 


bi 


280. 


Br(6i -> bxt) = 40% 


LI 


b 2 


354. 


Br(b 2 -> bxi) = 20% 



Table I: Masses of sbottoms, the lightest neutralino, and branching ratios of 61,2 — ► bxi f° r the test points considered in our 
numerical simulations. 



III. EVENT SIMULATION AND TEST POINTS 

We considered three scenarios for the new particle spectrum. Our first test point is the LHC point 5 (S5) [l(| that 
exhibits rather heavy sbottoms with sizeable decays into bxi; see Table I. Our second reference point is SPSla [ll| 
where the sbottom masses are close to ones in the first scenario, however the decays &i,2 — * bXx are suppressed. The 
third test point exhibits a somewhat light squark spectrum which could be eventually produced at a 1 TeV International 
e + e~ Linear Collider; we denote this point by LI. The masses for the LI parameter point are = 280 GeV, 
m^ 2 = 354 GeV, = 339 GeV, = 371 GeV, and m^o = 97.8 GeV. For the SPSla and S5 parameter choices 

the lighter of the two sbottoms is almost entirely a left state b\ ~ b^ while for LI it is approximately a maximal 
mixture of left and right states b\ ~ ^/2/2bL + \/2/2&r. The salient features of these reference points is summarized 
in Table I. 

In order to correctly treat all spin correlations we performed a parton-level analysis including the UED interactions 



in MADGRAPH [12j and usingSMADGRAPH [13j for the SUSY simulation. In our calculations we used CTEQ6L1 
parton distribution functions [14j with the factorization and renormalization scales are fixed by /j,f — Hb, — { m l + 
m^ 2 )/2 for the signal and /xp = /i/?, = my + X^Pt; for the SM backgrounds, where my is the mass of the relevant 
electroweak gauge boson and J2i Pr» 1S * ne sum of the transverse momentum of additional jets. 

We simulated experimental resolutions by smearing the energies (but not directions) of all final state partons 
with a Gaussian error. We considered a jet resolution AE/E = 0.5/a/B 0.03 and ajs T = 0.46 ^J2 e t for the 
missing transverse energy where ^2 Et is the sum of the jet transverse energies. In addition to that we included a 
b-tagging efficiency of e = 60% with a mistagging probability of 1/200 for light jets No K-factors were applied 
to signals or backgrounds since we do not expect QCD corrections to change significantly the shape of kinematical 
distributions [2, [3j. Nevertheless, given the increasing interest in studying spins correlations at the LHC by means of 
long and short decay chains, it is important to check if this is indeed the case. 



IV. ANALYSIS AND RESULTS 

At the LHC we analyzed the production of sbottom pairs followed by their decay in a b and the LSP 

PP 61,261,2 bbxlxi -> btyr , (8) 

i.e. the signal is characterized by two 6-tagged jets and missing transverse energy. One feature of the reaction {8j) is 
that the s-channel subprocesses qq — > bb* present the well-known angular distribution 

-^-ocl-cos 2 r (9) 

where 9* is the polar angle of produced scalar particles in their center-of-mass frame. However, this clean signature 
of the sbottom spin is contaminated by the subprocess gg — ► bb* which contains i-channel diagrams and quartic 
couplings; see Fig. [T] left panel. Therefore, it is easier to decipher the new state spin if we enhance the importance of 
the qq s-channel subprocesses via a judicious choice of cuts. Notwithstanding, the cuts to isolate the signal must not 
introduce bias in the angular distributions used to study the sbottom spin. 
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Figure 1: Left panel: cos 8* distribution for the production of sbottoms coming from qq and gg fusions. Right panel: same 
distribution as in the right panel but for the production of We used the test point S5 spectrum. 



On the other hand, the center-of-mass angular distribution of KK bottoms in UED produced by qq fusion is 

" Vr ocl+i ^ 1 "J } cos 2 9* , (10) 



dcosd* 




where M bl and E bl are the mass and energy respectively of the b^' in the center-of-mass frame. This distribution 
peaks in the forward and backward directions being quite distinct of the SUSY prediction. Moreover, we must also 
include the gg — > which contains t- and s-channel contributions which present a peak towards the forward and 

backward regions as well; see Fig. [2] right panel. 

At the LHC we can not reconstruct the the polar angle (#*) of produced particles in their center-of-mass frame 
due to the presence of undetected \\ or 7W. Therefore, we must use an alternative variable that retains part of the 
information carried by 9* . A convenient variable to use in our analysis is [f| 

cos%=tanh(^») (11) 

where Arj bb is the rapidity separation of the 6-tagged jets. Notice that A^f, is invariant under boosts along the 
collision axis, and consequently, cos 9 bb is invariant under boosts as well. The angle 9 bb is the polar angle between 
each reconstructed bottom jet direction in the longitudinally boosted frame in which the rapidities of the bottoms are 
equal and opposite. 

In order to cos 9 bb carry information of the produced particle spin, the flight directions of sbottoms and bottoms 
must be correlated. We depict in the left panel of Fig. [5] the cosine of the opening angle between the bottom and the 
sbottom flight directions in the bb* center-of-mass system for the S5 and LI spectra. Clearly, the bulk of the signal 
is concentrated in region of small opening angles as a consequence of the large energy of the sbottoms after cuts. 

Fig. [3] shows that cos 9 bb is indeed strongly correlated to the cosine of production polar angle (9*) of the sbottoms 
and £>W 'g in their center-of-mass system. Therefore, we must expect the shape of cos 9 bb distributions to resemble the 
polar angle spectra of sbottoms and KK bottoms apart from some smearing effects due to non-perfect correlations 
between the bottoms and sbottoms (KK bottoms) flight directions. Nevertheless, a clear separation between SUSY 
and UED distributions should be possible as was demonstrated in the case of smuon pair production at the LHC [f|. 
Taking a closer look at Fig. [3] we already realize that UED events are slightly more concentrated near cos6^ b = ±1 
while SUSY events are homogeneously distributed along the direction cos0£ 6 = cos 6^. 

In our analysis, we included the following backgrounds for the process pp — > bbfix'- 

• SM QCD and electroweak production of bbZ with Z — > vv that accounts for ~ 91% of the total background 
after cuts for the S5 and SPSla test points and ~ 72% for the LI scenario. 
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Figure 2: Cosine of the opening angle between the bottom and the sbottom in center-of-mass system for two different mass 
spectra considered here. 
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Figure 3: In the left (right) panel we display a scattered plot cos 6l b ( 
parameter point S5. 



i9* for sbottoms (b^ 1 ') production. Here we chose the 



• Reducible electroweak and QCD backgrounds like jjZ, jjW, jrv Tl bbj, jj, jjj, and jjb where some of the 
decay products evade detection and j stands for a light jet that might be mistagged as a b jet. These dangerous 
backgrounds are efficiently reduced by the requirement of two 6-tagged jets. The second largest background 
after cuts is bbW where the lepton from the W boson has a large rapidity \rji\ > 2.5. 

• SUSY processes, excluding bb* production, that lead to the final state bbxiXi- 

In order to properly trigger [l(| the event and tag the b-jets [l5[ the following acceptance cuts were imposed in all 
cases 



\Vb\ < 2.5 



p b T > 100 GeV 



AR bb > 0.7 



i> T > 100 



(12) 



A potentially large background is the QCD production of dijets once we take into account that mismeasurements 
of the jets properties can lead to missing transverse momentum and that this process has a huge cross section. This 
background can be efficiently reduced by the missing transverse momentum cut Eq. (fT2j) and by requiring that the 
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azimuthal angle between the jets and the missing transverse momentum satisfy 

\A$y T ,p Tj )\>0.3 . (13) 

In the LI test point simulations we applied the following cuts not only to enhance the signal and deplete the 
background but also to augment the importance of the qq s-channel subprocess 

il%>3G0GeV , M cff > 600 GeV , \A$(£r,p Tjsoft )\ < 2.4 , (14) 

where M e s is the sum of all jet and missing transverse momenta and A^>(p , T,PTj ao f t ) the azimuthal angle between 
the missing transverse momentum and the softest jet. After cuts the SUSY signal cross section is 38.1 fb while the 
background cross section is 2.4 fb, leading to S/B ~ 16.6. 

For the parameter point S5 we used the following cuts to optimize the signal 

M b5 > 600 GeV , M efi > 1 TeV , \A$ty T ,pT jaoft )\ < 2.4 . (15) 

The SUSY signal cross section at this test point is 4.55 fb with a background of 1.24 fb and S/B ~ 3.7 while the 
qq fusion accounts for ~ 40% of the signal events. On the other hand, ~ 45% of the UED signal stems from qq 
fusion. The cuts employed in this case were harder for two reasons: first, the bottoms are harder in this case once 
the sbottom is much heavier than the neutralino compared to LI case; second, the signal cross section is significantly 
smaller than in the LI case which required a deeper suppression of backgrounds to avoid a severe bias on the signal 
angular distributions. 

Finally, we imposed the following cuts for the test point SPSla 

M bl > 600 GeV , M eS > 1 TeV , (16) 

which lead to a signal (background) cross section of 1.07 (1.46) fb and S/B ~ 0.73. 

In the left panel of Fig. [H we show the impact of cuts on the da /d cos bb distributions for SUSY and UED predictions 
assuming the S5 spectrum. We normalized the UED signal cross section to the SUSY one. As the background events 
populate the bins of larger | cos 6 bb | , the selection cuts tend to suppress that region enhancing the central region 
for both SUSY and UED. The variable A<&(^x,PTj) is efficient in rejecting dangerous SM backgrounds like dijet 
production, however, it has a potential to bias the distributions as we can see in the left panel of Fig. [4] Therefore, a 
harder cut on this variable is not recommended in the study of spin correlations. 

A natural question at this point is whether we can mimic the SUSY signal by varying the UED mixing angles. 
According to Section [TT] there is not very much room to modify the UED Lagrangian to bring kinematical correlations 
closer to the SUSY prediction. The KK weak mixing angle 0w in Eq. ([5]) is fixed by the KK masses, so we can 
not change it while keeping the masses fixed. The same limitations hold when we try to adjust the mixing between 
the singlet and doublet KK fermions, described by the angle a^; see Eq. (|2|). In contrast to the 3rd-generation 
sfermion sector in the MSSM, the UED mixing angle is not a (third) free parameter, even if we move around the 
masses invoking boundary conditions. Notwithstanding, for illustration purpose we vary in the right panel Fig. 0] 
to check whether the SUSY cos can be reproduced by a UED decay chain with different couplings to the fermions. 
From Eq. ([6]) we see that varying effectively enhances the left or right couplings of the KK bottom decay into 
bottom plus LKP. This figure allows us to see that the changes in the UED parameters are not sufficient to mimic 
the SUSY predictions. 

We depict in Fig. [5] the cos 9^ spectrum with/without adding the background for the S5 test point, an integrated 
luminosity of 300 fb , and after applying cuts (|12|) and (|15[) . To avoid using any information but the spin we assume 
the S5 spectrum for the UED particles and normalize their production cross section times branching fractions to the 
SUSY rate. From Figure we can easily see that the production of fermionic strongly interacting states (UED) 
favor large separations (cos 6 bb ) between the ^-tagged jets while the production of scalar particles (SUSY) leads to 
a rather constant distribution. Note that these distributions indeed resembles the distributions of the production 
angles in the center-of-mass system which reflects the correlation between these observables. Moreover, the UED 
cos 9 bb distribution is similar to the SM background one since they take place through similar diagrams containing 
KK partners or SM particles with the same spin quantum numbers. 

For the assumed integrated luminosity it is rather easy to distinguish between the two models. These two possibilities 
can be disentangled, for instance, through the asymmetry: 

<7(|cosg; 6 | <0.5)-q(|cos^| > 0.5) 

cr(] cos0fc|< 0.5) +a{\ cos 6* bb \ > 0.5) ' 1 ' 
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Figure 4: Left panel: impact of cuts on the da / 'dcos8 bb distributions for SUSY and UED signals. We normalized the UED 
signal cross sections to the SUSY ones using the S5 spectrum. Right panel: SUSY and UED do jd cos bb distributions without 



backgrounds but varying the KK bottom mixing angle a' 1 '. 
M bb > 600 GeV. 



The basic cuts are the acceptance denned in Eq. (|12fl plus 




cosB , . cos® , 



Figure 5: In the left (right) panel we plot the cos9 bb distribution without (with) the addition of the background contribution. 
Here we used the parameter point S5 and assumed an integrated luminosity of 300 fb -1 . The error bars indicate the expected 
statistical uncertainties. 



This asymmetry is +0.238 ± 0.023 for the UED spin assignment while the SUSY interpretation it is significantly 
larger +0.373 ± 0.022 where the quoted errors are statistical. We estimate that an integrated luminosity of 300 fb -1 
is needed to reach a 5<r level signal for the S5 spectrum. 

The determination of the sbottom spin is much easier for the test point LI due to the large b — > bx± branching 
ratio and production cross section. We depicted in Fig. [6] the cos 6^ distribution and without/with adding the SM 
background. In this case a mere integrated luminosity of 15 fb -1 is enough to discriminate at the 5a level UED 
and SUSY. The asymmetries are given by +0.565 ± 0.019 for SUSY and +0.365 ± 0.021 for UED for this integrated 
luminosity. Therefore, there is a clear distinction between the two spin assignments (UED x SUSY) for both S5 and 
LI mass spectra and the above integrated luminosities. 

The results for the reference point SPSla are quite different from the S5 and LI ones. Due to small b — * bxl 
branching ratio, the SM backgrounds play a major role. We can see in the left panel of Fig. [7] that the pure SUSY 
and UED cos6^ b spectrum are quite different. However, once we add the SM background, which has a shape similar 
to UED, it is no longer easy to separate the SUSY from UED, even for an integrated luminosity of 1 ab _1 ; see the 
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Figure 6: Same as Fig. [5] but for the parameter point LI and an integrated luminosity of 50 fb . 
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Figure 7: Same as Fig. [5] but for the parameter point SPSla and an integrated luminosity of 1 ab . 



right panel of Fig. [7] For instance, the asymmetry (JTJI) is 0.282 ± 0.019 for SUSY and 0.200 ± 0.019 for UED and this 
extremely large luminosity. In this case SUSY and UED can be discriminated only at at ~ 4<r level. Therefore, it is 
important for this point to subtract the SM background which can be estimated, for instance, from the measurement 
of the bb/j, + fj," cross section. If we neglect systematic errors associated to this extraction from data, the sbottom 
spin can be determined at 5a level for an integrated luminosity of ~ 500 fb -1 . The left panel of Fig. [7] displays the 
distributions of SUSY and UED after a naive background subtraction where we just did not add the background to 
the signals and the error bars calculated as y/S + B. A dedicated study is necessary to determine the actual impact 
of the statistical and systematic errors [l7l |. 

In general, we noticed that the expected shape of distributions in cos are similar to the ones in smuon pair 
production and decay to LiXi @> which shows the universality and robustness of the method. Moreover, in both cases 
(b*sbx and the S5 spectrum seems to be more promising as compared to the SPSla spectrum, for example, 

SUSY versus UED discrimination of smuon spin assignments is possible in the case of S5 (SPSla) for an integrated 
luminosity of 200 fb -1 (500 fb" 1 ). 

As a final remark note that, apart from negligible effects from bb initiated contributions which include interactions 
with gluinos and electroweak interactions with Z bosons, the production rate does not depend upon the left or right 
nature of the sbottom since the QCD interactions to the gluons are blind to these details. On the other hand, the 
sbottom-bottom-neutralino vertex is sensitive to the relative content of mass eigenstates in terms of left and right 
states which by its turn depend upon the mass spectrum. As we have pointed out in the Section IIIII the lightest 
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sbottom is almost entirely a left state in the S5 mass spectrum, while it is an equal mixture of left and right states 
in the LI mass spectrum. As can be seen in the Figs. [5] and [H it seems plausible to conclude that the distributions 
are not sensitive to the particular mixtures of left and right states of the sbottom mass eigenstates, however a more 
detailed study including more test points is necessary to confirm this indication. 

V. CONCLUSIONS 

In the near future the LHC will start its endeavor in searching for new physics signals. Once those signals have been 
identified as the production of new states the next logical step will be the determination of the underlying model among 
all candidates by measuring the properties and interactions of the new particles. The size of production cross sections 
and mass spectra might provide valuable hints about the underlying model, nevertheless, an undisputed discrimination 
will only be possible after determination of the spins of the new states. We showed that the discrimination between 
a SUSY spin interpretation against an UED one is possible in the case of scalar bottom (fermionic KK bottom) pah- 
production in the reaction pp — > bbfix- Using the variable proposed in Ref. see Eq. we demonstrated that a 
clear determination of the spin of the decaying strongly interacting particle is possible provided the production cross 
sections and branching ratios into bxl are sufficiently large. Even in the worst scenario studied, the SPSla spectrum, 
the determination of spins might be possible after background subtraction, however, a dedicated study of the impact 
of statistical and systematic errors on the spin determinations will be necessary in this test point. 
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